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Comment on: 'Universal Behavior of Load Dis- 
tribution in Scale-free Networks' 

In a previous Letter [1], Goh et al have presented a- 
numerical study of the load — or betweenness centrality — 
distribution in a scale-free network whose degree distri- 
bution follows a power law p{k) ^ k^'^ where 7 g [2, oo[ 
is a tunable parameter. They showed that the load ^ is 
distributed according to a power- law P{e) ~ with 
exponent 5. On the basis of their numerical results, they 
conjectured that the value of i5 ~ 2.2 is independent of 
7 for the interval [2,3). Based on this apparent univer- 
sality, a classification of scale-free networks according to 
the value of i5 ~ 2.2 (class I) or 5 = 2 (class II) was pro- 
posed [2] . In this comment we argue that the value of 5 is 
not universal and varies significantly as 7 changes in the 
interval [2,3). The power law fits of the cumulative func- 
tion Prob(load > t) for the model proposed in [1] gives 
the values 5 = 1.84 ±0.04, 2.05 ±0.05 and 2.25 ±0.05 for 
7 = 2, 2.5, and 3 respectively, while for the BA model [3] 
5 — 2.3±0.1. The variations of (5 are significant enough to 
claim that it is not universal but in order to double-check 
our results we use an indirect way of computing 5. We 
study the relation between the load and the connectivity 
[1,4] which is of the form i ^ where the exponent 77 
depends on the network. As can be seen on Fig. (la), the 
power law holds remarkably for a large range of k allow- 
ing for an accurate measure of rj. We also checked that 
the value of 77 does not change significantly for different 
values of the system size (For 7 = 2.5, we obtain a rela- 
tive variation due to size going from N = 10^ to 5.10^ less 
than 1%). The exponents 77 and S are not independent 
and it is easy to show that [4] 77 = (7 — 1)/ (5 — 1). If the 
value of (5 ~ 2.2 is universal then 77 is a linear function 
of 7 with slope ~ 1/1.2 0.83. In Fig. (lb) we plot the 
measured 77 versus 7 for the different types of networks 
studied and the corresponding value predicted by univer- 
sality. This Fig. (lb) shows that if for 7 ~ 3 the value 
(5 = 2.2 seems to be acceptable, the claim of universality 
for 7 G [2,3) proposed in [1] does not hold (our results 
do not fit in the other class S = 2.0 either). In addition, 
we tested the universality for different values of 771 and 
we also obtain variations ruling it out: For 7 = 2.5 and 
for N = 2.10'^, we obtain 77 = 1.477 ± 0.006, 1.56 ± 0.006, 
and 1.64 ± 0.01 for 771 — 2, 4, 6 respectively. 
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FIG. 1. (a) Log-Log plot of the normalized average load 
versus connectivity for the same models as in [1] with m = 2. 
The power law fits (straight lines) give 77 = 1.27 ± 0.01 
{N = 3.10''), 1.467 ± 0.006 (iV = S.W), and 1.68 ± 0.02 
(iV = 5.10") for 7 = 2, 2.5, and 3 respectively. For the BA 
model, T] = 1.81 ± 0.02 (iV = 5.10*). (b) rj versus 7. If the 
universality proposed in [1] would be correct, the measured 
values for 7 £ [2, 3) should lie on the "universal" straight line 
corresponding to 5 = 2.2 (class I). 

The important exponent thus appears to be rj and it is 
interesting to note that 77 is significantly smaller than 
the maximum value rj — 2. This maximum value is 
reached when nodes with large centrality (ie. with large 
k) link together disconnected parts of roughly the same 
size. The load of these nodes is then of the order of 
I ^ k{k - l)/2 ~ k"^. The fact that 7; < 2 indicates that 
the different parts are also connected by shortest paths 
which do not pass through the central node. More gener- 
ally, it would be interesting to understand how 77 depends 
on the different parameters of the network such as 7 and 
the degree correlation. 

Marc Barthelemy 

CEA-Service de Physique de la Matiere Condensee 
BP12 Bruyeres-Le-Chatel, France 
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